Flux-creep in the second magnetization peak of BaFe1.9Nio.1As2 superconductor. 
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Flux-creep data was obtained for fields along the second magnetization peak observed in M(H) 
curves of BaFei.gNio.i As2 for if||c-axis. J/||a6-planes and H forming a 45° angle with afr-planes. The 
M-H loops from the different field directions can be collapsed onto a universal curve with a scaling 
factor equivalent to the superconducting anisotropy, showing that the pinning is three dimensional, 
although with remarkable differences in the vortex-dynamics as a function of field orientation. The 
resulting relaxation rate, 7?, when plotted as a function of field and temperature does not show any 
specific feature in the vicinity of the second magnetization peak field Hp, the relaxation shows a 
maximum at a field H2 well above Hp for H\\c and a minimum at H* for fields well below Hp for 
H\\ab and H-45°-ab. Isofield plots of the scaled activation energy obtained from flux-creep data at 
several different temperatures also do not show any evidence of a change in the pinning mechanism 
as the peak field is crossed. The Hp lines in the resulting phase diagrams do not appear to be 
consistent with a description-terms of a collective-plastic pinning crossover. 

PACS numbers: 74.70.Xa,74.25.Uv,74.25.Wx,74.25.Sv 



The study of vortex-dynamics in the novel pnic- 
tide superconductors 1 has attracted increasing attention, 
due to the considerable high-T c of these compounds, 
when compared to the conventional superconductors, 
and because of similarities with the high-T c cuprates 
superconductors 2 . Most of the pnictides systems exhibit 
the second magnetization peak, or fish-tail, in isothermic 
magnetization curves, as well as large flux-creep, allowing 
the study of different regions of the vortex-phase diagram 
in certain detail 3 - - — . The second magnetization peak is 
associated with a maximum in the critical current when 
measured as a function of field at a fixed temperature. 
This phenomena is not yet completely understood in 
pnictides^. Although pnictides have similarities with the 
cuprates, such as the layered structure and antiferromag- 
netism of the precursor non-superconducting system 2 -, it 
is well established that superconductivity in pnictides has 
a multi-band characther— for which it is predicted the ex- 
istence of non-Abrisosov vortice s 12 ' 13 which might lead 
to new effects in the vortex matter. In this work we 
study the vortex dynamics in a BaFe1.9Nio.1As2 super- 
conductor single crystal, by means of isofield magnetiza- 
tion M(H) curves and magnetic relaxation M(t) curves. 

Magnetization data were obtained by using commercial 
magnetometers: a 5T MPMS based on a superconducting 
quantum interference device (SQUID) was used for most 
of measurements with i/||ao-planes ; and a 9T PPMS was 
used for the other measurements including all data with 
-ff||c-axis and and i?-45°-a6-planes. The measurements 
were made after lowering the sample temperature from 
above T c in zero applied magnetic field (ZFC-procedure). 
The studied sample, is a high-quality single crystal of 
BaFe1.gNio.1As2 with transition temperature T c = 20 K, 
transition width AT C =0.3 K, mass = 120 mg and dimen- 
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FIG. 1: a) double plot of M(H) at 15 K for H\\c and H\\ab.; 
b) M(H) curves at 16 K for the three geometries. Insets: a) 
and upper b) shows a reduced plot of the same correspondent 
M{H) curves.; lower b), J C (H) at 14 K for H\\c. 



sions 0.5x1.6x0.02 cm. Details of the sample preparation 
can be found in Ref. [H]. For fz~||c-axis and H-A5°-ab ge- 
ometries, we carefully broke the sample and used a 43.1 
mg piece with dimensions 0.55x0.5x0.02 cm. For ff| lab- 
planes, the sample was attached to a hard-plastic slab 
perfectly inserted along the entire length of the straw 
tube used in the measurement systems, assuring an al- 
most perfect alignment of field with the a6-planes. 

Figure la shows a double plot of M(H) curves obtained 
at 15 K for H\\c and H\\ab where it is possible to observe 
the differences in the shape of the fish-tail between both 
field directions, which as shown below, are due to differ- 
ences in the vortex dynamics. Figure lb shows M(H) 
curves at 16 K for all field directions, all exhibiting the 
second magnetization peak, where Hon represents the 
onset of the second magnetization peak with its maxi- 
mum value occurring at Hp. The inset of Fig. la and 
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the upper inset of Fig. lb show a reduced plot of the 
curves appearing in the respective main figures, where 
M is divided by the respective value at Hp and H by 
the respective irreversible field Hirr. These inset figures 
show that Hirr and Hp, for the different geometries, can 
be collapsed onto a universal M-H curve. The lower in- 
set of Fig. lb shows a curve of the critical current J C (H) 
at T = 14 K for H\\c, which was obtained from the corre- 
spondent M(H) curve by using the Bean Model^. The 
almost perfect symmetry of the M(H) curves of Fig. 1 
with respect to the x-axis evidence that bulk pinning is 
dominant. Also, the equilibrium magnetization, Meq, 
defined as the average value of M on both branches of a 
M (H) curve, is very small, so we may use values of M 
instead (M — Meq) in the analysis that follows. Near T c 
for temperatures above 19 K, the second magnetization 
is no longer observed. 

The vortex dynamics study was performed by collect- 
ing isofield magnetic relaxation data, M(t) curves, along 
several isothermic M(H) curves, for fields lying below 
and above the second magnetization peak, for the three 
geometries. We also obtained isofield magnetic relax- 
ation curves as a function of temperature. Magnetic re- 
laxation data were collect for 2 hrs when obtained in 
both branches of M(H) curves and for 3.5 hrs when 
only in the lower branch. We also measured long time 
magnetic relaxation for approximately 12 hrs in differ- 
ent regions of the M(H) curve for ff||c-axis. All M(t) 
vs. log(i) curves show strictly linear behavior, starting 
above a transient time To ~ 3-4 min. Such a large tran- 
sient time was observed before in BaFei.82Nio.isAs2 with 
T c =8 K, and seems to be intrinsic of the material. The 
linear behavior with log(i) was also observed for the 12 
hrs relaxation curves. This fact suggests that it is more 
appropriate to analyse the data by using the relaxation 
rate R = dM/dln(t) as defined in Ref. 16. As shown 
in Ref. [IH one may obtain information on the apparent 
activation energy, but this quantity may not have any 
physical meaning in our data, since the magnetization 
Mq at the time t — 0, above which logarithmic relaxation 
should start, is not well defined, due to the ks 4 min long 
transient region (see the inset of Fig. 2c below). 

Figure 2 shows results of the analysis of M(t) data 
for if||e-axis. The logarithmic relaxation of the magne- 
tization, is exemplified by the selected curves plotted in 
the inset of Fig. 2c. Figure 2a shows a double plot of 
R{H) vs. H and the corresponding M(H) curves where 
the 2hrs relaxation data are also plotted. In Fig. 2a, we 
plot R as a function of increasing and decreasing field and 
we define H\ as the minimum in R (i? increases above 
Hi) and Hi a maximum (R decreases above H 2 ). It is in- 
teresting to note that the resulting R(H) curve resembles 
the correspondent M(H) curve after a shift to the right. 
The same trend as shown in Fig. 2a was observed for T 
— 14 K and 16 K data. To explore whether the field Hi 
is related to Hon, we took 3.5 hrs relaxation data in the 
lower branch of M{H) curves at T = 6, 8, and 10 K in- 
cluding more data for fields below Hon. Figure 2b shows 
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FIG. 2: Jf||c-axis: Double plots of —R vs. H and M(t) vs. H 
at a) T = 12 K; b) T = 8 K; c) R vs. H for the lower branch 
of M(t,H) curves. Upper inset: Double plot of R vs. H and 
M(t, H) for YBaCuO at T = 60 K. The lower inset shows 
selected magnetic relaxation data. 



the results for T = 8 K with a double plot of R vs. H 
and the correspondent M(t, H) vs. H , evidencing that 
the field Hi is indeed related to the field H on, but H 2 is 
about 30 kOe above Hp and these two fields are not likely 
to be related. It is clear from Figs. 2a and 2b, which are 
representative data from all experiments that there is no 
apparent change in the relaxation rate as the peak field 
Hp is crossed. Figure 2c shows the results of R vs. H for 
relaxation data obtained on the lower branch of all M (H) 
curves, where the arrows pointing up and down show the 
approximately positions of H2 and Hi respectively. For 
the sake of comparison, the upper inset of Fig. 2c shows 
a double plot, as in Fig. 2b, for a YBaCuO sample (T c 
~ 92 K), with relaxation data obtained during 60 min 
(if 1 1 c- axis) at the lower branch of an M(H) curve at T 
= 60 K. It is clear the perfect matching between Hp and 
the field position of the maximum in R, which in the 
case of YBaCuO 17 , represents a pinning crossover taking 
place as Hp is crossed. A direct comparison of the upper 
inset of Fig. 2c with Fig. 2a or 2b, suggests that the sec- 
ond magnetization peak in these two systems arise from 
different mechanisms. 

Figure 3 shows the results for if ||a6-planes. In that 
case, R decreases as field increases above Hon, reaching 
a minimum at a field denominated H* , which lies well 
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FIG. 3: R vs. H for the lower branch of M(t, H) curves for 
H ||afe-planes. Insets: Double plots of R vs. H and M(t,H)\ 
upper, at T = 16 K with H\\ab; lower, at T = 14 K with 
H45°ab: 



below Hp. For each curve of Fig. 3, arrows pointing up 
shows the position of Hp and pointing down of H*. A 
similar result was previously observed for an overdoped 
BaFei.82Nio.igAs2 with T c = 8 K for if ||c-axis and H\\ab- 
planes&. The upper inset of Fig. 3 shows a double plot of 
R vs. if and the correspondent M(t, if) vs. if curve. To 
check whether this change in R only occur when if | lab- 
planes, we measure the sample for an intermediate geom- 
etry, with H forming a 45° with a6-planes. The results 
for T — 14 K is shown in the lower inset of Fig. 3 (the 
same trend is observed for T = 16 K), showing that vor- 
tex dynamics for this geometry is similar to that observed 
for if 1 1 ab-planes. 

To check whether a pinning crossover, or a vortex 
phase transition, would perhaps become evident near Hp 
using a different approach, we obtained a few isofield 
magnetic relaxation data as a function of temperature 
for if||c and ff||ab. Now, as in RefUl, we calculate 
the pinning activation energy for a set of isofield data 
M(t, T) using the expression, U(M)=-Tln(dM/dt)+CT 
where C=hi(Bwa/ ird) is a constant, where B ps if is the 
magnetic induction, oj is an attempt frequency, a is the 
hop distance and d is the sample size. It is believed 18 
that the isofield U(M,T)/g(T/T e ) (where g(T/T c ) is 
an appropriated scaling function of U—) should be a 
smooth function of \M\ within a temperature region 
with same pinning mechanism. Figure 4 shows the re- 
sults of U(M,T)/g(T/T c ) plotted against \M\ for H\\c, 
as obtained for if = 25 kOe data with C = 10 and 
g(T/T c ) = (1 - T/T c ) 3 / 2 . Similar values of the con- 
stant C were found for cuprates superconductors^ and 
for pnictides^£. The almost perfect log|M | fit linking all 
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FIG. 4: U(M,T)/(1-T/Tc) 3/2 vs. |M| for H = 25 kOe, if He- 
axis. The solid line represents a log\M\ curve. Inset: R vs. T 
for H = 20 and 25 kOe, H||c-axis. 



the data in Fig. 4 suggests the existence of only one pin- 
ning mechanism over the entire AT range, for which Hp 
is located near 11 K. The same trend was observed for 
U(M,T)/(1 - T/T c ) 3 / 2 vs. \M\ obtained for if \\ab data 
with if = 15 kOe (not shown). The inset of Fig. 4a shows 
plots of R vs. T as obtained from isofield data with if = 
20 and 25 kOe for if||c, which do not show any visible 
effect near ifp located at 12 K and 11 K in each respec- 
tive curve. One would expect some feature in the plot of 
R vs. T as Hp is crossed, either for a pinning crossover or 
for a vortex-lattice phase transition. The same behavior 
was observed on similar plots for if ||a6-planes with if = 
15 and 20 kOe. 

Figure 5 shows the resulting phase diagram where 
characteristic fields of the three geometries are 
plotted against temperature after being divided 
by ^f{sin{9)) 2 + (l/3)(cos(6>)) 2 where 6 is the an- 
gle between if and the afr-plane, and the factor 
3^Hp(H\\ab)/Hp{H\\c)KiHirr(H\\ab)/Hirr(H\\c) is 
of the order of the system anisotropy. 19 For clarity 
reasons we did not plot values of ifi which lie close to 
Hon. The collapse of the if on, Hirr and Hp lines is 
evident and each line follows a (1 — T/T c ) m dependence 
with m = 1.8, 1.6, and 1.4 respectively. Figure 5 
also suggests a collapse of if* values ( for if \\ab and 
H45°ab) which lie between Hon and Hp. Values of 
if2 (for if||c) lie between Hp and Hirr. One may 
associate the anisotropic vortex dynamics observed 
here with the anisotropic neutron spin resonance found 
on a similar optimally doped sample^, which is cor- 
roborated by isotropic vortex dynamics 8 - and isotropic 
neutron spin resonance^ found on overdoped samples of 
BaFe2_ :E Ni l; As2. Although we could not find evidence 
for a pinning crossover associated with Hp, we notice 
that the peak field position Hp is time dependent, as 
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FIG. 5: Phase diagram for the three geometries after scale the 
Y-axis for H\\ab and H45°ab. Solid lines represents a fitting 
of the data. Dashed lines are only a guide to the eyes. 

predicted within the plastic pinning scenario^. For this 
reason, we try to fit the Hp line for each geometry by the 
proper expression^ Hpm(l — (T/T c ) 4 ) 14 which failed. 
This fact concur with the experimental evidence that, 
Hp in the studied system, is not related to a pinning 
crossover. 

Although our results do not show a direct evidence for 



a vortex lattice phase transition, neither the literature 
seems to present studies of the vortex structure for in- 
termediate and high fields for the studied system, the 
strong evidence of absence of a pinning crossover associ- 
ated with the peak effect motivated us to explore the pos- 
sibility that the peak effect in the critical current might 
be due to a softening of the vortex lattice associated with 
a phase transition as predicted in Ref.22. The lower in- 
set of Fig. Ib shows a fitting of J C (H) at 14 K for H\\c 
to the expression J e (B)=A/[(B - Bp) 2 + (AB) 2 } 5 / 4 of 
RefH3 where A is a fitting parameter Bp is the peak po- 
sition and AB is the peak width. The resulting fitting 
conducted in a wide field range is excellent and produced 
Bp = 14 kOe and AB = 24 kOe. Similar fittings were 
obtained at different temperatures with H\\c. 

In conclusion, our study of flux-creep data based on 
the rate of magnetic relaxation, and on isofield activation 
energy curves, did not show any evidence of a pinning 
crossover occurring near the second magnetization peak 
of M(H) curves in BaFe1.gNio.1As2. Curves of R vs. H 
for H \\ab and i/T5°a&-planes show a minimum at H* well 
below Hp, which behavior is quite similar to that found 
for an overdoped crystal of the same system for H\\ab and 
H ||c£. Curves of R vs. H for H \\c show a somewhat differ- 
ent behavior with a minimum at H\ associated to Hon, 
and a maximum at H2 well above Hp. A fitting based 
on the expression predicted by plastic pinning failed to 
explain the Hp line of each phase diagram. 
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